Properties of high- Tc copper oxides from band models of spin-phonon couphng. 
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The mechanism of spin-phonon couphng (SPC) and possible consequencies for the properties of 
high- Tc copper oxides are presented. The results are based on ab-initio LMTO band calculations and 
a nearly free-electron (NFE) model of the band near Ef- Many observed properties are compatible 
with SPC, as for the relation between doping and q for spin excitations and their energy dependence. 
The main pseudogap is caused by SPC and waves along [1,0,0], but it is suggested that secondary 
waves, generated along [1,1,0], contribute to a 'waterfall' structure. Conditions for optimal Tc, and 
the possiblities for spin enhancement at the surface are discussed. 

PACS numbers: 74.25. Jb,74.20.-z,74.20.Mn,74.72,-h 
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I. INTRODUCTION. 

The normal state properties of high-Tc copper ox- 
ides show many unusual properties like pseudogaps, 
stripe-like charge/spin modulations with particular en- 
ergy and doping dependencies, Fermi-surface (FS) "arcs" 
in the diagonal direction, 'kinks' and 'waterfalls' (WF) 
in the band dispersions, manifestations of isotope shifts, 
phonon softening and so on Band results for 

long 'l-dimensional' (1-D) supercells, calculated by the 
Linear Muffin-Tin Orbital (LMTO) method and the lo- 
cal spin-density approximation (LDA), show large spin- 
phonon coupling (SPC) within the CuO plane of these 
systems This means that an antiferromagnetic 

(AFM) wave of the correct wave length and the proper 
phase is stronger when it coexists with the phonon 
The LMTO results have been used to parametrize the 
strength of potential modulations coming from phonon 
distortions and spin waves of different length ^13i] . These 
parameters have been used in a nearly free electron 
(NFE) model in order to visualize the band effects from 
the potential modulations in 2-D. Many properties are 
consistent with SPC, as have been shown previously 
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II. CALCULATIONS AND RESULTS. 

Ab-initio LMTO band calculations based on the 
local density approximation (LDA) are made for 
La(2_a;)Baa:Cu04 (LBCO), with the use of the virtual 
crystal approximation (VCA) to La-sites to account for 
doping, X. Calculations for long supercells, mostly ori- 
ented along the CuO bond direction, are used for mod- 
eling of phonon distortions and spin waves [l^ . The cal- 
culations show that pseudogaps (a 'dip' in the density- 
of-states, DOS) appear at different energies depending 
on the wave lengths of the phonon/spin waves. This is 
consistent with a correlation between doping and wave 
length, and with phonon softening in doped systems p^ . 

The difficulty with ab-initio calculations is that very 
large unit cells are needed for realistic 2D-waves. An- 
other shortcoming is that the original Brillouin zone is 



folded by the use of supercells, which makes interpreta- 
tion difficult. The band at Ep is free-electron like, with 
an effective mass near one, and the potential modulation 
and SPC can be studied within the nearly free-electron 
model (NFE) J^]. The AFM spin arrangement on near- 
neighbor (NN) Cu along [1,0,0] corresponds to a potential 
perturbation, V{x) = V*exp{—iQ ■ x) (and correspond- 
ingly for y). A further modulation (q) leads to ID-stripes 
perpendicular to x (or "checkerboards" in 2-D along x 
and y), with a modification; V{x) = V*exp{~i{Q — q)-x), 
and the gap moves from the zone boundary to {Q — q)/2. 

The NFE model reproduce the qualitative results of 
the full band (1-D) calculation. In 2-D it leads to a cor- 
relation between doping and the amplitude of V^, be- 
cause the gap (at Ep) opens along {kx,0) and {0,ky), 
but not along the diagonal [l3]. The combined effect 
is that the dip in the total DOS (at Ep) will not ap- 
pear at the same band filling for a small and a wide gap. 
The q vs. X behavior for a spherical NFE band with m* 
close to 1, and parameters (for one type of phonon) 
from ref. [l3| show a saturation, see Fig.l. This is quite 
similar to what is observed Q. The reason is that no 
checkerboard solutions are found for larger doping, but 
unequal qxfi (fixed near 0.11) and qy/2 produce real- 
istic solutions. The DOS at Ep, N, is lowest within 
the pseudogap. A gap caused by spin waves, disappears 
at a temperature T* when thermal excitations can over- 
come the gap, and the spin wave can no longer be sup- 
ported. Therefore, the spin-waves are most important 
for the pseudogap (even though phonons are important 
via SPC), and T* is estimated to be l/4th of the spin 
part of [ijj. The opposite a;- variations of T* and N 
(note that A oc NV^) provides an argument for opti- 
mal conditions for superconductivity for intermediate x 
p^ . Moreover, the pseudogap competes with supercon- 
ductivity at underdoping, since dynamic SPC would be 
their common cause. However, there is a possibility to 
raise N and Tc by creating an artificial, static pseudogap 
through a periodic distribution of dopants or strain. Two 
parameters, the periodicity and the strength of the per- 
turbing potential, should be adjusted so that the peak in 
the DOS above or below the static pseudogap coincides 
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FIG. 1: The variation of magnetic modulation vector g/2 (Cir- 
cles), density-of-states in the pseudogap (+-signs N{Ef)/50 
in states per cell/Ry/spin), and T*/10000 (x-signs in K), as 
function of doping a; in a NFE model with parameters from 
ref. J.4]. Two different g- vectors along qx and qy are needed 
for doping larger than about 0.12 as indicated by the filled 
and open circles. 
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FIG. 2: Filled circles: Calculated q — hw relation from the 
2D-NFE model and the parameters Vg for doping x =0.16. 
The solution without SPC at the largest energy, is less precise. 
Broken line: Approximate shape of the experimental disper- 
sion as it is read from figure 3c in the work of Vignolle et al 
[lO| for La2-xSr3,Cu04 at a; = 0.16. Thin semi-broken line: 
Experimental dispersion read from the data by Tranquada et 
al [y| on LBCO at lower doping, x = 0.125. 



with Ep- 

The degree of SPC is different for different phonons. 
The total with contributions from phonons and spin 
waves, calculated from LMTO and information from 
phonon calculations for Nd2Cu04 [16], are 17, 18, 23 and 
22 mRy at energies centered around 15 (La), 25 (Cu), 50 
(plane-0) and 60 meV (apical-0), respectively p^]. The 
results for a;=0.16 are shown in fig. 2 together with ex- 
perimental data d, [l3| • The points below 70 meV are for 
the coupling to the 4 types of phonons. The spectrum is 
shaped like an hour-glass with a "waist" at intermediate 
energy with largest SPC for plane-0. The solutions for 
energies larger than 70 meV are independent of phonons 
and the exact (5, w) behavior is more uncertain [1^ . Less 
doping implies larger V* and longer waves. All g become 
smaller and the waist becomes narrower, as can be ver- 
ified in LBCO for x = 1/8 0], and recently in lightly 
doped Lai.g6Sro.o4Cu04 '4'|. However, the spin modula- 
tion in the latter case is in the diagonal direction. Heavier 
0-isotopes will decrease the frequencies for the phonons 
and the coupled spin waves, and move the waist to lower 
E. 

Another odd feature is the WF-dispersion of the band 
below Ep in the diagonal direction, seen by ARPES [3. 
The suggestion here is that this feature comes from a gap 
below Ep in the diagonal direction. An inspection of the 
potential for stripe modulations along [1,0,0] reveals that 
the potential becomes modulated also along [1,1,0], albeit 
in a different fashion. The potential is slowly varying 
like the absolute value of sme-functions with different 
phase along different rows. No NN-AFM potential shifts 



are found along [1,1,0], and the dominant Cu-d lobes of 
wave function for k along [1,1,0] are oriented along [1,1,0] 
and not along the bond direction. Various arguments for 
the effective periodicity, partly based on these conditions, 
indicate that a gap should appear at about 1/3 of the 
distance between T and the M point when the doping 
is near 0.16. The effective Vg should be less than half 
of the amplitude along [1,0,0]. The result is shown in 
Fig. 3. The k-position of the gap and the extreme values 
of the gap energies 0.5-1 eV below Ep) are not too 
far from what is seen experimentally [3], but again, the 
quantitative power of the NFE-model is limited. It is not 
clear if the vertical part of the band dispersion can be 
observed. A vertical line is connecting the states above 
and below the gap in Fig. 3, which could be justified for 
an imperfect gap away from the zone boundary. 

The dynamics is important if SPC mediates supercon- 
ductivity. But static, stripe like features are identified 
by surface tunneling spectroscopy (STM) [13]. Impuri- 
ties and defects near the surface might be important, 
but also the surface itself could modify the conditions for 
SPC. The latter hypothesis is investigated in LMTO cal- 
culations which simulate the surface through insertion 
of two layers of empty spheres between the outermost 
LaO-layers. These calculations consider 3 and 5 layers of 
undoped La2Cu04, and 3 layers of a doped LBCO with 
and without phonon distortion in a cell of length 4ao in 
the CuO bond direction. The SPC remains in the surface 
layer. The effective doping is in all cases found to increase 
close to the surface, which has 0.1-0.2 more electrons/Cu 
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Ep. In addition, also the apical-0 nearest to the surface 
acquires a sizable moment. This calculation is simplified, 
with a probable interaction across the empty layer, but 
it shows that static AFM surface configurations are very 
close to stability. 
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FIG. 3: Thin lines: LMTO band structure for LBCO between 
M and T. Broken line: the FE fit, and the heavy hne the NFE 
solution with a gap. Ep is at zero in undoped LBCO, and at 
the thin broken line for x ~ 0.15. 

than the Cu in the interior, and the magnetic moment 
is 2-3 times larger in the surface layer. The moments 
disappear without field, but a calculation for 3 layers of 
La2Cu04 with a narrower separating layer, has stable 
AFM moments ±0.06/xb per Cu within the surface layer, 
and the local DOS on the Cu at the surface drops near 



III. CONCLUSION 

Band calculations show that SPC is important for 
waves along [1,0,0] or [0,1,0], with secondary effects in 
the diagonal direction. Many properties, like pseudo- 
gaps, phonon softening, dynamic stripes, correlation be- 
tween q and x, smearing of the non-diagonal part of the 
FS, and abrupt disappearance of the spin fluctuations at 
a certain T*, are possible con seq uences of SPC within a 
rather conventional band Different SPC for dif- 

ferent phonons leads to a hour-glass shape of the {q, lu)- 
spectrum with the narrowest part for the modes with 
strongest coupling. The much discussed WF-structure in 
the diagonal band dispersion could be a result of a sec- 
ondary potential modulation in this direction and a gap 
below Ep. Static potential modulations within the CuO- 
planes, such as for superstructures, could compensate the 
pseudogap and enhance N{Ep) and Tc- Spin waves be- 
come softer through interaction with phonons and near 
the surface. These LDA results show a tendency for static 
spin waves at the surface. 
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